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Abstract 

The Nil geometry, which is one of the eight 3 -dimensional Thurston 
geometries, can be derived from W. Heisenberg's famous real matrix group. 

The aim of this paper to study lattice coverings in Nil space. We intro- 
duce the notion of the density of considered coverings and give upper and 
lower estimations to it, moreover we formulate a conjecture for the ball ar- 
rangement of the least dense lattice-like geodesic ball covering and give its 
covering density A 1.42900615. 

The homogeneous 3-spaces have a unified interpretation in the projec- 
tive 3-sphere and in our work we will use this projective model of the Nil 
geometry. 



1 Notions of the Nil geometry 

In this Section we summarize the significant notions and denotations of the Nil 
geometry (see [[I], [|6l). 

*Mathematics Subject Classification 2010: 52C17, 52C22, 52B15, 53A35, 51M20. 
Key words and phrases: Thurston geometries, Nil geometry, lattice-like geodesic ball coverings. 
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The Nil geometry is a homogeneous 3-space derived from the famous real 
matrix group L(]R) discovered by Werner Heisenberg. The Lie theory with the 
methode of the projective geometry makes possible to investigate and to describe 
this topic. 

The left (row-column) multiplication of Heisenberg matrices 



1 X z\ / 1 a 
1 y\ io 1 b 
.0 1/ \0 1. 



'1 a + X c + xb + 
1 b + y 

,0 1 



(1.1) 



defines "translations" L( 



= {{x,y,z) : X, y, z e R} on the points of the 
space Nil = {(a, b, c) : a, b, c e M}. These translations are not commutative in 
general. The matrices K.{z) < L of the form 



'1 z^ 
K(z) 3 \ 1 I ^ (0,0,^) 
,0 1, 



(1.2) 



constitute the one parametric centre, i.e. each of its elements commutes with all 
elements of L. The elements of K are called fibre translations. Nil geometry 
of the Heisenberg group can be projectively (affinely) interpreted by the "right 
translations" on points as the matrix formula 



1^1; a, b, c) (1; a, b, c) 



(I X y z\ 

10 

1 X 

\^o ly 



{l]x + a,y -\-b,z -\-bx -\- c) (1.3) 



shows (see (1.1)). Here we consider L as projective coUineation group with right 
actions in homogeneous coordinates. We will use the Cartesian homogeneous 
coordinate simplex EQ{eo)E'^{ei)E^{e2) E^{e^), ({ej} C with the unit 
point £"(6 = eo + ei + 62 + 63)) which is distinguished by an origin Eq and by the 
ideal points of coordinate axes, respectively. Moreover, y — ex. with < c e IR 
(or c G M \ {0}) defines a point (x) = (y) of the projective 3-sphere VS^ (or 
that of the projective space where opposite rays (x) and (— x) are identified). 
The dual system {(e*)}, ({e*} C V^) describes the simplex planes, especially 
the plane at infinity (e°) = Ef'E^Ef, and generally, v — defines a plane 
(w) = (v) of VS^ (or that of V^). Thus = xw = yv defines the incidence of 
point (x) = (y) and plane (w) = {v), as (x)l(i<,) also denotes it. Thus Nil can be 
visualized in the affine 3-space (so in E^) as well. 
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Tlie translation group L defined by formula (1.3) can be extended to a larger 
group G of coUineations, preserving the fibering, that will be equivalent to the 
(orientation preserving) isometry group of Nil. 

In [2] E. Molnar has shown that a rotation trough angle cu about the z-axis at 
the origin, as isometry of Nil, keeping invariant the Riemann metric everywhere, 
will be a quadratic mapping in x, y to 2;-image z as follows: 



7^ = r(0, u) : (1; x, y, z) (1; x, y, z); 
X = X cos u — ysinu, y = xsinu + y cos u, 

: = z — -xy + -{x'^ — y"^) sin 2uj + -xy cos 2u. 



(1.4) 



This rotation formula 71, however, is conjugate by the quadratic mapping M. to 
the linear rotation i7 in (1.5) as follows 

M : (1; X, y, z) ^ (1; x' , y', z) = (1; x,y,z- ^xy) to 



n-. il;x',y',z')^il;x\y\zn = il; x', y', z') 



with : (l;x",i/",z") 



M- 



/l 0\ 

cosw sinw 

— sin cos uj 

\0 1/ 

1 



(l;x,y,^) = (l;x",y",z" + -x"y"). 



(1.5) 

This quadratic conjugacy modifies the Nil translations in (1.1), as well. Now a 
translation with (X, Y, Z) in (1.3) instead of [x, y, z) will be changed by the above 
conjugacy to the translation 



{l]x,y,z) — > {l]x,y,z) = {l;x,y,z) 



(\ X Y Z-\XY\ 

2 
2 ^ 



1 
1 
\0 



(1.6) 



that is again an affine coUineation. We shall use the following important classifi- 
cation theorem. 

Theorem 1.1 (E. Molnar ||2l) (1) Any group of Nil isometries, containing a 3- 
dimensional translation lattice, is conjugate by the quadratic mapping in (1.5) to 
an affine group of the affine (or Euclidean) space = E"^ whose projection onto 
the (x,y) plane is an isometry group ofE?. Such an affine group preserves a plane 
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— )■ point polarity of signature (0, 0, ±0, +). 

(2) Of course, the involutive line reflection about the y axis 

(l;x,y,z) (l;-x,y,-z), 

preserving the Riemann metric, and its conjugates by the above isometrics in fl) 
(those of the identity component) are also Nil-is ometries. Orientation reversing 
Nil-isometry does not exist. 

The geodesic curves of the Nil geometry are generally defined as having lo- 
cally minimal arc length between their any two (near enough) points. The equation 
systems of the parametrized geodesic curves g{x{t),y{t), z(t)) in our model can 
be determined by the general theory of Riemann geometry: We can assume, that 
the starting point of a geodesic curve is the origin because we can transform a 
curve into an arbitrary starting point by translation (1.1); 

x{0) — y{0) — z{0) — 0; x{0) — coos a, y(0) = csina, 
i(0) = w; — TT < a < TT. 

The arc length parameter s is introduced by 

s = Vc^ + • t, where w = sin6', c = cos6', — ^ < 
i.e. unit velocity can be assumed. 

Remark 1.1 Thus we have harmonized the scales along the coordinate axes. 

The equation systems of a helix-like geodesic curves g{x{t),y{t), z{t)) if < 
\w\ < 1: 

2c . wt /wt \ , , 2c , wt , /wt 



, , /c wt /wt \ ic wt (Wt \ 

x[t) = — sm — cos I — + a ) , y{t) = — sm — sm — + a 1 , 

r r/. sin(2wi + 2q;) — sin2Q;\ 

Zit) ^ wt ■ \l + —: 1 ^ — ^ + 

^' I 2^2 LV 2wt J 

/ sin(2wt)\ / sin(wi + 2q;) — sin2Q;\ 1 1 

'^t / V 2^t JJj 

( c^ r/^ sin(u't)\ /l-cos(2u't)\ . , ^ 



(1.7) 



In the cases w = the geodesic curve is the following: 

x{t) = c ■ tcosa, y{t) = c ■ tsina, z{t) = -^(? ■ t^ cos asm. a. (1.8) 
The cases |w| = 1 are trivial: (x, = (0, 0), z — w - t. 
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Definition 1.1 The distance d{Pi, P2) between the points Pi and P2 is defined by 
the arc length of the geodesic curve from Pi to P2. 

1.1 On the geodesic ball 

In our work [|6l we have introduced tlie followin definitions: 

Definition 1.2 The geodesic sphere of radius R with centre at the point Pi is 
defined as the set of all points P2 in the space with the condition d{Pi,P2) = 
R. Moreover, we require that the geodesic sphere is a simply connected surface 
without selfintersection in the Nil space. 

Remark 1.2 We will see that this last condition depends on radius R. 

Definition 1.3 The body of the geodesic sphere of centre Pi and of radius R in 
the Nil space is called geodesic ball, denoted by Bp^{R), i.e. Q G Bp^{R) iff 
< d{Pi,Q) < R. 

Remark 1.3 Henceforth, typically we choose the origin as centre of the sphere 
and its ball, by the homogeneity o/Nil. 

We apply the quadratic mapping M : Nil — > at (1.5) to the geodesic 
sphere S, its A^-image is denoted by S' = Ai{S). 

We choose a point P{x{R, 6, a),y{R, 9, a), z{R, 6, a)) lying on a sphere S 
of radius R with centre at the origin. The coordinates of P are given by param- 
eters (a G [-vr,7r), 6 G [-f,f], R > 0) (see (1.5), (1.10)), its yW-image is 
P'{x'{R, e, a), y'{R, 9, a), z'{R, 9, a)) G S' where 



+ a 



2c . 


wR 


/wR 


= — sm 




[ 2 


w 


2 




2c 


wR . 1 


rwR 


= — sm 


2 '^^( 




w 




V 2 



y'{R, 9,a) = — sin sin — h a 

c^R c^ -K TT (1-9) 

z'{R, 9,a)=wR+^-£-^^ sinwR, {9e[--,-]\ {0}), 

if 6* = then x'{R, 0, a) = Rcosa, 
y'{R, 0, a) = Rsin a, z'{R, 0, a) = 0. 

We can see from the last equations that {x'Y + {y')"^ = ^ sin^ ^ and that the 
^'-coordinate does not depend on the parameter a, therefore S' can be generated 
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by rotating the following curve about the z axis (lying in the plane [x, 2;]): 

TA/n n\ 2c wR 2cos6' RsinO 
X{R, e) = —sm—= . sm — - — , 
w 2 smd 2 

c^R 

Z(R, 6)=wR + -smwR = 

^ ' ' 2w 2w^ (1.10) 

+ Tita^ - ^"'(''^"''')- <" 1-2- 2I \ 

if ^ = then 0) = R, Z{R, 0) = 0. 



Remark 1.4 From the definition of the quadratic mapping Ai at (1.5) it follows 
that the cross section of the spheres S and S' with the plane [x, z], is the same 
curve which is specified by the parametric equations (1.10). 

Remark 1.5 The parametric equations of the geodesic sphere of radius R can be 
generated from (1.10) by Nil rotation (see (1.4)). 

We have denoted by B{S) the body of the Nil sphere S and by B{S') the 
body of the sphere S', furthermore we have denoted their volumes by Vol{B{S)) 
and Vol{B{S')), respectively. 

In H we have proved the the following theorem: 

Theorem 1.2 The geodesic sphere and ball of radius R exists in the Nil space if 
and only if R & [0, 27r]. 

We obtain the volume of the geodesic ball of radius R by the following integral 
(see 1.10): 



— 

VoliBiS))=2ix 



2 cos 9 . (R sin ^) \ ^ / 1 -R cos'^ 9 cos 9 sin {R sin 1 



= 271 I sin ^ • ^ h 

io V sin^ 2 / V 2 sin^^ sin^ 

cos^ 9 sin (R sin 9) IR cos^ 9 cos [R sin 9)\ 

^ ^9 2 ^9 ) 

(1.11) 

The Nil sphere of radius R is generated by the Nil rotation about the axis z 
(see the equation system (1.10) and remarks (1.4), (1.5)). The parametric equation 
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system of the geodesic spliere S{R) in our model: 

, 2 cos e 



z{R,e,(j)) = wR + 
RsmO + 



2c . 


wR 


= — sm 




w 


2 


2c 


wR 


= — sm 




w 


2 



sin 9 

2 cos 9 



Rsm9 
sm ■ cos ( 



2 

i?sin^ 



sin 9 



sm ■ 



sm I 



c^_ 

2w 
cos^ 9 



cos (, 
■ sin (, 

C2 

— - sin wR H — ( — sin | sm : 

2^2 A\w 2 

• /r^ • n\ 1 /2cos0 . „sin6'\2 
sm(itsmt^) + - ( — : — ^smit 



l/2c . wi?\2 



2 sin 6* 2 sin^ 61 ' 4 V sin 6* 

— IT < ^ TT, - 

if ^ = then x{R,0,(f)) = R cos 0, y{R,0,(j)) = R sin , 
z{R,0, 



-<9<-2ind9^0. 
2 - - 2 ^ 



sin 20 



(1.12) 



1 2 

- R cos sin ( 



The following theorem was obtained by the derivatives of these parametrically 
represented functions (by intensive and careful computations with Maple through 
the second fundamental form) (see [l6l): 

Theorem 1.3 The geodesic Nil ball B{S{R)) is convex in affine-Euclidean sense 
in our model if and only if R E [0, |]. 



1.2 The discrete translation group L(Z, k) 

We consider the Nil translations defined in (1.1) and (1.3) and choose two arbi- 
trary translations 



(I t\ t\ tl\ 

10 

I t\ 

\0 1 / 



and T2 



(\ ^2 ^2 '^iX 

10 

1 

\0 1 / 



(1.13) 



now with upper indices for coordinate variables. We define the translation (rs)'^, {k 
G N, A; > 1) by the following commutator: 



/I 

1 

1 

\0 






1 



(1.14) 
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If we take integers as coefficients, their set is denoted by Z, then we will generate 
the discrete group ((ri, T2), k) denoted by L(ti, T2, k) or by L(Z, k). 

We know that the orbit space Nil/L(Z, k) is a compact manifold, i.e. a Nil 
space form. 

Definition 1.4 The Nil point lattice Tp{ti, T2, k) is a discrete orbit of point P 
in the Nil space under the group L(ti, T2, k)= 'L{'L, k) with an arbitrary starting 
point P for all (/c e N, /c > 1). 

Remark 1.6 For simplicity, we have chosen the origin as starting point, by the 
homogeneity o/Nil. 

Remark 1.7 We can assume that t{ — 0, i.e. the image of the origin by the 
translation Ti lies on the plane [x,z]. 

In the following we investigate the most important case k — 1 where T3 correspond 
to the fibre translation, i.e. = t2^Ti^T2Ti. 

We illustrate the action of L(Z, 1) on the Nil space in Fig. 1. We consider 
a non-convex polyhedron — OT1T2TST12T21T23T213T13, in Euclidean sense, 
which is determined by translations ti, T2, T3. This polyhedron determines a solid 
T in the Nil space whose images under L(Z, 1) fill the Nil space just once, 
i.e. without gap and overlap. 

Analogously to the Euclidean integer lattice and parallelepiped, the solid T 
can be called Nil parallelepiped. 

T is di fundamental domain of L(Z, 1). The homogeneous coordinates of the 
vertices of T can be determined in our affine model by the translations (1.13) and 

(1.14) with the parameters t\, ^ = 1 (i e {1, 2}, j e {1, 2, 3}) (see Fig. 1 and 

(1.15) ). 

ri(i,ti,o,t?), T2{\AA^% T3(i,o,o,^), 

Ti3(l, t\, 0, ^ + t\), Ti2(l, t\ + 4, t\, tl + 

1 1 2 1 2 3 3 1 2 3 At\ ^^'^^^ 
221(1, ^1 + ^2) ^2) ^1^2 + ^1 + ^2)) ^23(1, ^2) ^2) ^2 "I 

T213 = T23i(l, t\ + t\,tl, (A; + 1)^ + + tl). 
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Figure 1: 



We liave determined in [|6l the volume of by the following integral: 

Voli^) = \t\-tl\ dxdy = it\ ■ tlf. (1.16) 

Jo 

From this formula it can be seen that the volume of the Nil parallelepiped depends 
on two parameters, i.e. on its projection onto the [x, y] plane. 

2 The lattice-like geodesic ball coverings 

A family of subsets /C = {Ki)i(zi of Nil, / is a set of indices is called covering 
of Nil if each point of Nil belongs at least one of the set {Ki, i G /}, i.e. 
Nil = [Ji^j Ki. A covering of Nil space is a lattice covering if it is of the form 
K. = (Kq + y)veL{z,k) where Kq is a element of {Ki, i G 1} i.e. the lattice 
coverings those coverings which cover Nil by translated copies of a single body 
Kq and in addition the translates are vectors of a lattice L(Z). 

In following, we are only considering lattice coverings consisting of geodesic 
balls of the Nil. Let (R) denote a geodesic ball covering of Nil space with balls 
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B'^{R) of radius R where their centres give rise to a Nil point lattice r(Ti, T2, 1). 
J^o an arbitrary Nil parallelepiped of this lattice (see (1.13), (1.14)). The images 
of J'o by our discrete translation group L(ri, T2, 1) cover the Nil space without 
overlap. 

Remark 2.1 In the d-dimensional Euclidean space Yfi, {d > 1) an arbitrary 
lattice r(ri, T2, 1) under the group L(Z, 1) gives a lattice covering of equal balls 
if the radius R of the balls is large enough, but this is not true in the Nil space, 
because of a geodesic ball exists in the Nil space if and only if R e [0, 27r] (see 
Theorem 1.2). 

If we start with a lattice covering Bf^{R) and shrink the balls until they finally do 
not cover the space any more, then the threshold value of the shrinking radius R 
defines the least dense covering of equal balls to a given lattice r(ri, r2, 1). The 
thresfold value i?"^ is called covering radius of the point lattice r(ri, r2, 1): 

R'' := mm{R : where B^(R) lattice covering by r{ri, T2, 1)}. (2.1) 

For the density of the packing it is sufficient to relate the volume of the "min- 
imal covering ball" to that of the solid J^o- 

Analogously to the Euclidean case it can be defined the density 5{R'^, ri , T2, 1) 
of the lattice-like geodesic ball covering Br{R'')'. 

Definition 2.2 



The main problem is that to which lattice L(ti, T2, 1) belongs the minimal density 
Aopf We introduce for the "optimal arrangement" the following denotations: 



A{R',n,T2,l) : 



Vol{Br{R')) 
Vol{^) 



(2.2) 



^opt{R'opt,r^,rll) min | 



Vol{BriR''j) 
Vol{r) 



(2.3) 



Remark 2.3 The covering radius is the radius of the circumsphere of the lattice's 
Dirichlet-Voronoi poly tope, that is the largest distance between the midpoint and 
the vertices of its Dirichlet-Voronoi poly tope. 
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2.1 The lattice-like ball covering of the lattice Lopt{ri^\ T2^^ 1) 

First we consider tlie ball arrangement SJ?'* (Ropt) of the densest lattice-like geodesic 
ball packing in the Nil space, given by formulas (2.4), (2.5) (see [16]). 



(2.4) 



(a) d{0,T,) = 2R = d{T,,T,), 

(b) d{0,T2) = 2R = d{T2,T^), 

(c) rf(Ti,T2) =2i?, 
{d) d{0,T^) = 2R. 

^ 1.30633820, = Ropt, Ropt ~ 0, 73894461; 
4'°^* ^ 0, 65316910, tl'°P^ ^ 1, 13132206, ts'"^* ~ 1.10841692, (2.5) 

rpOpt /-| ,l,Opt n ,3,Opt\ rpOpt /i ,l,Opt ,2,Opt ,3,Opt\ 

This packing can be generated by the translations Ij optijl^^ , , 1) where r°^* and 
Tg^* are given by the coordinates i = 1, 2; j = 1, 2, 3 (see (2.5)). Thus we 
obtain the neigbouring balls around an arbitrary ball of the packing B'p'^{Ropt) by 
the lattice r(r°^*, t^^, 1). We have ball "columns" in 2;-direction and in regular 
hexagonal projection onto the [x, ?/] -plane. 

The Nil parallelepiped ^ = 0T°'"T^'t;^'T^'T^'T°1'T^IT^' is a fun- 
damental domain of L(Z, 1). The homogeneous coordinates of its vertices are 
known (see Fig. 1 and (1.15)). 

We examine the covering radius R^ to the lattice r(r°^*, t^^^ , 1). 

Rp := min{_R : where B^{R) lattice covering by r(r°^*, Tg^*, 1)}. 

It is sufficient to investigate such ball arrangements B^opt (R) where the balls 
cover ^ or the Nil solid V = OT°'''T°l^'T°'''TfTll'T°l'T°l' (see Section 1.2). 
From (2.4) and (2.5) follows, that the point sets {0, Tf^*, T^*, Tg"^*}, {T°'p\ 

rpOpt rpOpt rpOpt^ J rpOpt rpOpt rpOpt rpOpt^ J rpOpt rpOpt rpOpt rpOpt^ J rpOpt rpOpt rpOpt 

T°|*}, {T°2*, 7^21'*; 72f T""!*} are congruent by Nil isometrics. The radius R of 
each circumscribed ball to the above point sets can be determined by the following 
system of equation: 

d{0, C) = d{C, Tf) = d{C, T^') = d{C, T^*), (2.6) 

where C(l, c^, c^, c^) is the center of the circumscribed ball of the point set {0, T°*'*, 
T^*, Tg^*} (d is the Nil distance, see Definition 1.1): 

^ 0.45981062, ^ 0.26547179, ^ 0.79997799, 

(2.7) 

R ^ 0.90293941. 
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Remark 2.4 C{l,c^,c^,c^) is a vertex of the Dirichlet-Voronoi domain of the 
point O in the Nil space. 

R e [0, |] thus by Theorem 1.3 the ball of radius R is convex in affin-Euclidean 
sense. We form tetrahedra oT^^'T^^'T^'', T^^'TfT^'T^', t^^^t^'T^'t!^'", 

rpOptj,Optj.Optj,Opt^ joptrpoptr^optrpopt^ rp^trpOptr^OptrpOpt EUClidean SenSe, Which 

fill the Nil solid V = OT^^T^^T^^T^^T^^T^^T^* just once. Their circum- 
scribed congruent Nil balls are convex thus they cover the tetrahedra and so the 
ball arrangement B^opt{R) cover the Nil solid V. Thus the radius R of circum- 
scribed ball give us the covering radius R^ to the lattice Lopt(rf^*, Tj'^*) and we 
get by (1.11), (1.16) and by the Definition 2.2 the following results: 

Vol{B{Rl)) w 3.12538516, Vol{V) = Vol{F) w 2.18415656, 

A(i?^, rr , , 1) := M^I^M ^ 1.43093459. ^^'^^ 

It follows that 

^opt{Rl,t, rl 1) < A(i?^, rf*, rf*, 1) ^ 1.43093459. (2.9) 



Remark 2.5 The density of the least dense lattice-like ball covering in the the 
Euclidean space is 

/^opt{Rtpt,r^,r^, 1) < As = ^ « 1.46350307. 

2.2 Upper estimation for the covering radius 

We consider a arbitrary lattice covering B^{R'^) where F = r(ri, r2, 1) and 

R'^ := inf {i? : where B^iR) lattice covering by r(Ti, T2, 1)}. 

The fundamental domain T of the translations group Lopt(ri, r2, 1) and its ver- 
tices with their homogeneous coordinates can be seen in our affine model in (1 . 15) 
with parameters t^, i G {1, 2}, j e {1, 2, 3} (see Fig. 1). We divide the Nil solid 
V = OT{rx2T2T{rxzT2{r2z (have been derived from the fundamental domain f. 
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see Section 1.2) into Nil solids OT{r2T^, T^T{r2zT^z, T3T1T23T2, risTiTsaTa, 
T12T1T23T13, T12T21T23T13 which are tetrahedra in terms of Euclidean geometry. 
It is clear, that one of them contain the centre of its circumscribed Nil ball. Sup- 
pose now, that this "tetrahedron" is OT1T2T3 of which circumscribed ball B{R) 
centred by C passing through the points O, Ti, T2, T3. We note here that from 
conditions of the Nil ball follows that the ball B{R) contain the "Euclidean line 

segment" OT3 = \J Vol{T) (see Section 1.1-2). From the Definition 2.2 follows, 
that 

A{R ,Ti,T2,l) > A{R,Ti,T2,l) = Y^i(^^^ ^ OTI ■ ^ 

Remark 2.6 A geodesic curve in Nil space "parallel to the axis z" correspond 
to an "Euclidean line segment" (see (1.5), (1.7)) with same lenght. 

Farther estimation for the density we need to investigate the upper bound of 

the length of the line segment OT3 = \J~Vol{^ in B{R). 

Let (R) be a geodesic sphere of radius R with centre at the origin. We apply 
the quadratic mapping M. : Nil — > at (1.5) to the geodesic sphere S^, its M.- 
image is denoted by S'^ = Ai{S^), moreover we have denoted by -8(5'* (i?)) the 
body of the Nil sphere S^R) and by B{Si{R)) the body of the sphere Si{R), 

Lemma 2.7 The length of the vertical chords ofS^{R) do not change at M.. 

The proof of this lemma follows from the definition of the quadratic mapping M.. 



Lemma 2.8 B(S'^(R)) = J\4{B{S^{R))) is convex in our model in Euclidean 

sense if and only if R G [0, tt]. 

Proof: From the Section 1.1 can be seen that S'^{R) can be generated by rotating 
the curve €(6) = {X{R,e), Z{R,e)) {6 e [-f,f] \ {0}) (see 1.10) about the 
z axis (lying in the plane [x, z]) thus the convexity of the ball B{S'^{R)) follows 
from the investigation of the derivative 

dZ(R,e)) _ 

de ~ 

^ / (- (cos {e)fRsm (9) + 2 sin (sin (9) R) (sin {9)f) 

= cos(^)f . (2.11) 

^ 2(sm(t^)) 

((cos {9)f cos (sin {9) R) t sin (9) + 2 (cos {9)f sin (sin (9) R)) ■ 



2 {sin {9)y 
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From the first (2. 11) and second derivatives of Z{R, 9)) follows that if R E [0, vr] 
then the ball B{Sl{R)) convex. If R e {n,2n] then the equation ^^^^ = 
possesses a solution and the curve C has an inflection point in the interval 9 E 
(0, |). In Fig. 2 can be seen the complete curve C{9) = {X{R, 9), Z{R, 9)) for 
the parameter R = 27r, then the curve possesses at the point 6* = | the maximum, 
Z(27r,f) = f. 




Figure 2: 



Lemma 2.9 If Re [I, 2tx] then 



> A(i?;;, r°P*, rf *, 1) ^ 1.43093459. 



Proof: 



1. i?e [f,7r] 

We get by (2.10) and by Lemmas 2.6-7 the following inequalities: 



A(i? , r„ 1) > —^fT- > (2i?)2 • (2.13) 



The function /(-R) = ^"(2^)2^^'' depends only on the parameter R and 
its graph which is increasing on the interval R E [f,7r], can be seen in 
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Fig. 3. Consequently, at the point i? = | the function possesses a mini- 
mum, /(f) ^ 1.71179510 > A(i?;;,r°P*,rf*, 1) ^ 1.43093459. 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

1.75 2.0 2.25 2.5 2.75 3.0 



Figure 3: 



2. Re [tt, 27r] 

We divide this interval into two part, [tt, 2n] 



[7r,f]U(f,27r] 



(a) R G [vr, f ] 

In this interval the greatest vertical chord hi for the balls B(R) is 



h^ = l^r^ 10.21017613. 



(2.14) 



37r] 
' 2 i 



The Fig. 4. a show its increasing function on the interval R E [n 
It is clear that this function possesses a minimum, at the point R = tc. 

fi{n) ^ 1.441711246 > A(i?;;, *, r^*, 1) ^ 1.43093459. 

(b) i?e (f ,27r] 

Similarly to 2/a case we examine the possible greatest vertical chord 
/i2 for the balls B(R) on the given interval. We get that /i2 = 57r ^ 
15.70796327 (see Fig. 5). 

Vol{Br{R)) 



A{R%n,T2,l) > 



hi 



:= f2{R). 



(2.15) 
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The Fig. 4.b shows the graph of the f2{R) on the interval [^,27r]. 
It is evident that this function has a minimum, at the point R = 

/2(f) ^ 2.372757787 > A(i?;;, r^^*, rf *) ^ 1.43093459. 




Figure 4: a, b 

Immediate consequence of the Lemma 2.8 is the following Theorem: 

Theorem 2.10 The radius R^^^ which belongs to the minimal density Aopt {Ropt^ ^i^ 
r|, 1) of lattice -like geodesic ball coverings, is found in the interval [0, |]. 

Remark 2.11 The optimal covering Nil ball B{R^p^) is convex in Euclidean 
sense (see Theorem 1.3 and jj^). 

2.3 Lower bound to the covering density AoptiR^pt: rf , r|, 1) 

In this section we consider a arbitrary lattice covering Bf.{R'^) where F = F (ri , r2 , 1 ) 
(see (2.1-2)). The fundamental domain T of the translations group L(ri, r2, 1) is 
given by its vertices in our affine model in (1 . 15) (see Section 1.2 and Fig. 1). Sim- 
ilarly to Section 2.2 we divide the Nil solid P = OT1T12T2T3T13T21T23 (have 
been derived from the fundamental domain see Section 1.2) into Nil solids 

l-f 2-' 3? -'S-'l-' 23-^135 -'3-'1-'23-'2j 12-' l-^ 23-* 2j 12-' l-^ 23-* 13j 12-^ 21-^ 23-* 13 

which are tetrahedra in Euclidean sense. It is clear, that one of them contain the 
centre of its circumscribed Nil ball and thus can be assumed, that this "tetrahe- 
dron" is OT1T2T3 of which circumscribed ball B{R) centred by C(l, c\ c^, c^) 
and passing through the points O, Ti, T2, T3. By Theorem 2.10 can be supposed 
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that R G [0, |] thus B{R) is a convex geodesic Nil ball in Euclidean sense and 
contain the Euclidean tetrahedron OT1T2T3. 

In order to find a lower bound to the covering density we investigate the density 
function 

A{R,Ti,T2, 1) = _ , ~, — = (2.16) 



Vol{J^) OT^ 



to a given volume of parallelepiped Vol{J^) = OT^. We have to find the minimum 
radius of the circumscribed ball of the Nil solid OT1T2T3 if the Euclidean line 
segment OT-s is given. 

We project the points Ti and T2 parallel to the z axis onto the equidistance 
surface of O and T3 which is a hyperbolic paraboloid in our model with equation 
2z — xy = t\ ■ t^. Their images are Tf and T|. 

Lemma 2.12 The radius R^ of the circumscribed ball with centre C^(l, c^, c^, c^) 
of the Nil solid 0TlT2T^ which is a tetrahedron in terms of the Euclidean geom- 
etry, at most R. 

Proof: From the conditions of the Nil balls follows (see Section 1.1, [i6| and 
Fig. 2) that vertical line segments from the point upto the equidistance surface are 
contained by the circumscribed ball B(R), thus R^ < R. Moreover, during this 
projection the volume Vol{J') = OT^ = \t\ -tl^ does not change thus 

A{R,Ti,T2,l) > — = —2 . (2.17) 



Vol{J^) OT, 



Remark 2.13 Here we do not examine whether Bt{R^), F = r(rf , r|, 1) is a 
covering or not. 

We consider the Euclidean plane a which passes through the point and per- 
pendicular to the axis y. If we move in this plane then t\ is constant thus 
Vol{^) = OT^ = \t\ ■ t2|2 does not change. By Theorem 2.10 and by Remark 
2.1 1 we have obtained the following 

Lemma 2.14 The radius R^ of the circumscribed ball of the Nil solid OT1T2T3 
is minimal at the above moving of the point Tf if a touch the Nil ball B{R^) (see 
Fig. 5). 

Remark 2.15 At the "minimal position" of the point Tf is c^ = t^. 
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o'=t; 




X I 



Figure 5: 

Forthermore we shall decrease the radius of the circumcribed ball while the vol- 
ume of the parallelepiped is constant. We translate the Nil solid OTfTl^^T-i and 
its circumscribed ball with Nil translation further we apply the quadratic 

mapping M. : Nil — > at (1.5) to this arrangement (see Fig. 5). At these 
transformations the volume of the parallelepiped Vol{J') = OT^ and the vol- 
ume of the circumscibed ball Vol{B{RP)) do not change. Thus we get a solid 
0'T[T2T^ in affin space A'^ with its circumscribed ball of radius EF centred by 
origin. This ball is convex in Euclidean sense (see Theorem 2.10). Moreover 
the points T[, lie in the plane [x,y] and OT3 = 0% = lA{H'T{T^). Here 
we have denoted by H' E [x, y] the midpoint of the line segment O'T^ and by 
A{H'Tm) the area of the Euclidean triangle H'T[T^. 

Working in analogy with what we know from Euclidean geometry, if we fix 
the volume of the parallelepiped Vol{T) = OT^ = 0% = \ A{E'T[T'^ then by 
the Fig. 6. a,b,c,d can be derived the following 

Lemma 2.16 The radius EF of the circumscribed ball of the solid O'TITI^T^ to a 
given volume of parallelepiped is minimal if H'T[ = if'Tg where H'T[ and H'T2 
are Euclidean line segments. 

We have to examine these arrangements. To each RP can be determined the 
line segment O'Tg thus can be investigated the density function A(i?^) of these 
orders and we need to test this function for a possible minimum if E [0, |]. 
A(_RP) can be examined by careful computation with Maple. The graph of the 
function A(R^) can be seen in Fig. 7 and we get the following resuls: 

A(i?Ln) - 1.36278112, ^ 0.85847445. (2.18) 
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o'=t; 



o'=t; 



Ax , 















r' 


1 ► 





















o'=t: 




b. 



o'=t: 




Figure 6: 



Note tliat, it is easy to prove, that the ball arrangement belonging to the above 
given lattice does not yield a geodesic ball covering in the Nil space. 



Corollary 2.17 Consequently, we have obtained a lower bound to the covering 
density: 



AoAKpt^r^^r^, 1) > A(i?L„) ^ 1-36278112. 



(2.19) 
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Figure 7: a, b 

3 Conjecture for the least dense lattice-like 
ball covering in the Nil space 

First we consider a Nil lattice T which is generated by the translations L(ri, r2, 1) 
where ri and r2 are given by the coordinates t\ i = \, 2; j = 1, 2, 3 (see, Fig. 1). 
Moreover the points Ti and T2 lie on the equidistance surface of points O, T3 

and tl = tl = -^^-|^, and one immediate consequence, that \tl\ = 1^1-^21 = 

We have ball "columns" in 2;-direction and in regular hexagonal projection 
onto the [x, y] plane. 

Remark 3.1 The lattice Topt generated by the translations Lopt(r{'^*, r2^*, 1) (see 
(2.5)) is one of the above lattices. 

The radius R of the circumscribed ball of the Nil solid OT1T2T3 to a given pa- 
rameter t\ can be determined by the following system of equation: 

diO, C) = die, n) = d{C, Ti) = die, T2), (3.1) 

where C(l, c^, c^, c^) is the center of the circumscribed ball of the point set {0, Ti, 
T2,Ts}. 

In order to find the "suspected lower covering density " we investigate the 
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density function 



^ voKBrm ^ vomm 

Vol(T) \it\Y4V- 



To every R can be determined tlie parameter t\ tlius can be examined tlie density 
function A(i?) of these arragement and we need to test this function for a possible 
minimum if i? e [0, |]. A(i?) can be investigated by careful computation with 
Maple. The graph of the function A(i?) can be seen in Fig. 7.b and we get the 
following resuls: 

^(Rmin) ~ 1.42900615, ^ 0.86046718, (3.3) 



^i,min _ 1,26001585, ti'™" ^ 0.68746826; 
t^'™'^ ^ 0, 63000792, ^3'™" ~ 1, 09120574, t^'™" ^ 1.03120239, (3.4) 

rpmin /-I ,l,min n ,3,mm\ rp /-i ,l,mm ,2, mm ,3,mm\ 

Similarly to the Section 2.2 it is easy to see, that the ball arrangement belonging 
to the above given lattice is a geodesic ball covering in the Nil space, thus we get 
the following 



Theorem 3.2 



1.36278112 ^ A(Cm) < ^optiKpt, ^2^ 1) < A(i?mm) ^ 1.42900615. 



Conjecture 3.3 The least dense lattice-like ball covering in the Nil space is de- 
rived by the Nil lattice Tmin which is generated by the translations Lmini'r™"', 
r™", 1) where r™" and r™" are given by the coordinates tf"^^"' i = 1,2; j = 
1, 2, 3. The minimal covering radius R'^p^ = R'^in ~ 0.86046718 and 

Aopt{R'op„T^,T^,l) = A{R:^J ^ 1.42900615. 



Our projective method gives us a way of investigation the Nil space, which 
suits to study and solve similar problems (see [6 |). In this paper we have examined 
only some problems, but analogous questions in Nil geometry or, in general, in 
other homogeneous Thurston geometries are timely (see [Q, flUl, ifTOll ). 

Acknowledgement: I thank Prof. Emil Molnar for helpful comments to this 
paper. 



22 



Jend Szirmai 



References 

[1] Molnar, E. The projective interpretation of the eight 3 -dimensional ho- 

mogeneous geometries. Beitrage zur Algebra und Geometrie (Contri- 
butions to Algebra and Geometry), 38 (1997) No. 2, 261-288. 

[2] Molnar, E. On projective models of Thurston geometries, some rele- 

vant notes on Nil orbifolds and manifolds. Siberian Electronic Math- 
ematical Reports, \http://r^semr math.nsc. ru\ 7 (2010), 491-498. 



[3] Molnar, E. — Prok, L — Szirmai, J. Classification of tile-transitive 

3-simplex tilings and their realizations in homogeneous spaces. 
Non-Euclidean Geometries, Jdnos Bolyai Memorial Volume Ed. 
Prekopa, a. and Molnar, E. Mathematics and Its Applications 
581, Springer (2006), 321-363. 

[4] Scott, P. The geometries of 3-manifolds. Bull. London Math. Soc, 15 

(1983) 401-487. (Russian translation: Moscow "Mir" 1986.) 

[5] Szirmai, J. The optimal ball and horoball packings to the Coxeter hon- 

eycombs in the hyperbolic d-space. Beitrage zur Algebra und Geome- 
trie (Contributions to Algebra andGeometry), 48 No. 1, 35-47,(2007). 

[6] Szirmai, J. The densest geodesic ball packing by a type of Nil lattices. 

Beitrage zur Algebra und Geometrie ( Contributions to Algebra and 
Geometry), 48 No. 2, 383-398, (2007). 

[7] Szirmai, J. The densest translation ball packing by fundamental lat- 

tices in Sol space. Beitrage zur Algebra und Geometrie ( Contributions 
to Algebra and Geometry), 51(2) (2010), 353-373. 

[8] Szirmai, J. Geodesic ball packing in S^xR space for generalized Cox- 

eter space groups. Beitrage zur Algebra und Geometrie ( Contributions 
to Algebra and Geometry), to appear 2011. 

[9] Szirmai, J. Geodesic ball packing in H^xR space for generalized Cox- 

eter space groups. Mathematical Communications, to appear 2011. 



[10] Szirmai, J. Lattice-like translation ball packings in Nil space. 

Manuscript to Publicationes Mathematicae, Debrecen 2011. 



On lattice coverings of Nil space by congruent geodesic balls 



23 



[11] Thurston, W. P. (and Levy, S. editor) Three-Dimensional Geometry 

and Topology. Princeton University Press, Princeton, New Jersey, Vol 
1 (1997). 



